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Abstract: We present a general approach for coupling a specific mode 
in a planar photonic crystal (PC) waveguide to a desired free-space mode. 
We apply this approach to a Wl PC waveguide by introducing small 
index perturbations to selectively couple a particular transverse mode to 
an approximately Gaussian, slowly diverging free space mode. This "per- 
turbative photonic crystal waveguide coupler" (PPCWC) enables efficient 
interconversion between selectable propagating photonic crystal and free 
space modes with minimal design perturbations. 
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1. Introduction 

Photonic crystals (PCs) exhibit a photonic bandgap that enables localization and manipulation 
of light at the wavelength-scale. |jj. In two-dimensional (2D) photonic crystals, an important 
element is the photonic crystal waveguide, which enables compact optical delays (|2][3l, slow 
light and enhanced nonlinear processes BUS], and efficient connections between various on- 
chip devices iniiiiii. 

Various techniques exist for coupling to photonic crystal waveguides, including chip-to-fiber 
butt-coupling f9\, directional coupling from a tapered fiber ifTOlfTTI . and vertical grating cou- 
plers including slanted 1121 or focusing lens gratings lfT3l fT4l . Although these gratings can 
reach high efficiency, they constrain the PC device design as they disrupt the PC lattice and 
typically affect all waveguides modes. Here, we introduce a general approach for introducing 
highly non-invasive grating couplers that rely on small perturbations to existing device designs. 
For instance, in a Wl photonic crystal waveguide in a holey PC membrane, the grating consists 
only of small dielectric index perturbations to a set of holes near the waveguide region. We call 
this coupling device a Perturbative Photonic Crystal Waveguide Coupler (PPCWC). Although 
this approach is general, we focus here on the situation where we introduce extra holes adjacent 
to the existing air holes. This restricted set of perturbations eliminates fabrication challenges 
that may be associated with non-circular or very small, isolated perturbations. 

Similar directional coupling from photonic crystal cavities by introducing small perturbation 
has been explored in previous studies ITSl [161 [TtI . The index perturbation is derived here us- 
ing coupled mode theory, under the constraint that it can be easily integrated into the photonic 
crystal fabrication. This approach enables a straightforward design of perturbations that causes 
selective coupling between a desired PC waveguide mode and a desired free-space mode (cho- 
sen here to be a Gaussian). Figure[TJa) shows the coupling scheme. When light propagates along 
the waveguide (the red arrows) and passes through the coupler, the periodic index perturbations 
act as a grating and light scatters out of the slab. We use three-dimensional finite difference time 
domain (FDTD) simulations to calculate the waveguide modes for a given index perturbation 
and to find the efficiency of the resulting grating coupler We find that even very small per- 
turbations allow a high coupling efficiency between the free-space radiative mode and a given 
waveguide mode, while causing only very small changes to the propagation of other waveguide 
modes. As demonstrated in the end of Section 2, the directional scattering from this coupler is 
controllable by detuning the feed frequency or modulating the perturbation size. In Section 3, 
we discuss the mode selectivity of our device by investigating the problem in ^-space, and in 
Section 4 we analyze the transmittance and the modification of the shape and the directionality 
of the scattered beam. 



2. Coupled mode theory 



We consider a guided wave in a Wl photonic crystal line defect waveguide and the light scat- 
tered at the grating coupler We derive the coupled mode equations from wave equation 

where e(r) is the relative dielectric constant and c is the speed of light in vacuum. We define e„, 
as the relative dielectric constant for the waveguide, which is assumed to be periodic along the 
propagating direction. The solution of Eq. (1) with e = e„, is in form of 

E,, = B,,{r)e'^"''^^^'-^-^\ (2) 

B„ are Bloch states having periodicities as the multiples of the periodicity of £„, where 
k represents the crystal momentum and x is the propagating direction of the waveguide. We 
desire a scattered beam propagating in free space denoted by A(r). A is also a solution of Eq. 
(1) with the relative dielectric constant e = Eq. We assume the field to be a superposition of the 
waveguide and free space modes: 

E = a{t)A{r)e''^(^' + Y^Bi{r)e''^^'''^'[b{k,t)e-'^^^^^ (3) 

k 

Here, a{t) is the slowly varying component of the Gaussian beam, and b(k,t) and c(k,t) are 
the slowly varying components of the forward and backward propagating Bloch modes for the 
waveguide, respectively. Substituting Eq. (3) into Eq. (1), we can obtain the full coupled mode 
equations describing the interaction of the coefficients a{t), b(k,t), and c(k,t) through Ae. In 
particular, the coupling rate between a{t) and b(k,t) can be evaluated EKHl to be 

2fdr^P-\A\^ 

where ASn- = £i — £»., representing the index perturbations introduced into the already existing 
structures. The subscript letters w and t denote the existing waveguide and the target structure. 




Fig. 1. (a) Schematic showing the coupling system, (b) The dispersion relation of a Wl 
PC waveguide. The fundamental mode (the blue line) and a higher mode (the green line), 
having frequencies co = 0.2436 and 0.3 15 1 {2Ttc/a), respectively, at wave-vector kx = n/a, 
are labeled in the figure. 
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Fig. 2. FDTD Simulation of PC waveguide and calculated coupling strength of index per- 
turbed vertical coupler, (a) (b) x and y components of electric field of the waveguide mode 
CO = 0.2436. (c) A Gaussian beam profile which we expect as the scattered light, (d) Dot 
product of the y component of the waveguide mode and the target Gaussian mode, (e) 
Calculated coupling strength Ki,g (Eq. (4)) for target Gaussian modes with different diver- 
gences, (f) The resulting perturbed PC structure (the left panel) and the index perturbation 
distribution (the right panel) of the case with the perturbation size of 0.3244, in units of the 
fractional area of unit cell, in (e). 



From Eq. (4), the coupling constant is determined by A* ■ and the perturbation term Ae„,. 
Since the phase of A is nearly constant in the plane of the PC, whereas varies rapidly along 
the direction of propagation, the spatial integral of the dot product A* ■ B-^ will vanish. For a 
non-vanishing coupling constant, the index perturbation Ae,y should have the same periodicity 
as B^, assuming an interaction region that is large compared to the period a. 

Now we consider a realistic waveguide made by introducing a line defect in a hexagonal 
photonic crystal lattice with periodicity a, refractive index of the slab material n = 3.6, slab 
thickness h = 0.7 a, and hole radius r = Q.3a. The full coupling system is as shown in Fig.fTla). 
We used a 3D FDTD simulation to calculate the waveguide modes. Figure [Tfb) shows the dis- 
persion relation of the Wl photonic crystal waveguide. A fundamental and higher order mode 
inside the photonic bandgap (PBG), wedged by the areas shaded in gray, are shown in blue and 
green, respectively, and are confined below the light line, represented by the red line. Figures 
Wia) andlltb) show the x and y electric field components of the fundamental mode, which we 
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Fig. 3. Structures and simulation results of a line defect PC waveguides with and index 
perturbation introduced. Figures shows the distributions of the y component of the electric 
field of the waveguide modes at (a) CO = 0.2436 and (b) CO = 0.3078, respectively. 



label as Bj_, at wave-vector — n/a, and frequency (0 = 0.2436, in units of iKc/a, where c is 
the speed of light in free space. As shown previously, the index perturbation must have a period 
of 2n/kx = 2ax to scatter this fc-state to the F point in fc-space, where the in-plane = 0. Note 
that because the x and y components of the electric field of the waveguide mode have different 
symmetries along the y direction, an index perturbation with even symmetry across y = will 
scatter only the y-component, but not the .^-component of the field. Similarly, it is possible to 
scatter only the x-component with a anti-symmetric perturbation distribution across the y = 
line. For the following analysis, we will choose to scatter the y-component. Furthermore, we 
suppose that we want to scatter into a vertical mode that has a Gaussian profile in the in-plane 
directions, as shown in Fig. I2c). Figure |2d) shows the dot product A* ■ B-^ which appears 
in the integrand in Eq. (4). To find a perturbation that scatters light into a slowly diverging 
Gaussian beam, we calculate the coupling rate for beams with different divergence angles and 
compare the coupling constants among different index perturbation designs. We discuss in this 
paper, the perturbations are formed by enlarging three-hole clusters with periodicity 2a in the 
i-direction. The blue and green lines in Fig.|2e) represent the coupling constants Ki and K2 for 
Gaussian beams, focused in the plane of the photonic crystal, with waists of 2a and a, respec- 
tively, and the red line represents the ratio JCi/k^. The horizontal axis denotes the perturbation 
size in units of the area of the unit cell of the triangular lattice, and the vertical axes denote the 
coupling constant (left) and the ratio K1/K2 (right). We observe that a local maximum exists 
for Ki I K2 when the index perturbation size is 0.3244, in units of the fractional area of unit cell, 
that is, the light scattered in this case diverges less than a well-defined Gaussian beam. The 
perturbed structure for this case is shown in Fig. ^f) and the corresponding coupling strength 
Ki and the ratio ffi/fc^ are 0.2228 and 1.815, respectively. As expected, the coupling rate into 
the Gaussian beam of smaller waist is lower than for the wider, less rapidly diverging beam. 
The brown shaded area in the right panel of Fig. |2f) shows the index perturbation. 

We check our theoretical predictions using another 3D FDTD simulation of the coupling ef- 




Fig. 4. (a) The intensity distributions of thie y component of thie electric field, Ey, of the 
scattering from a 5-period long coupler, of which the perturbation pattern is as shown in 
|2lf), viewed in a vertical plane passing through the middle of the waveguide, where z = 
is the center of the PC slab. The waveguide is excited at a> = 0.2436. (b) The three 
dimensional view of the density isosurface of the y component of the electric field, Ey. The 
structure of the PC waveguide with a coupler is shown at the bottom. 

ficiency of a grating consisting of the index perturbation as shown in Fig.|2f). The fundamental 
waveguide mode at frequency ft) = 0.2436 is excited by a Gaussian beam source at one end of 
the PC waveguide. The y-component of the electric field in the plane of the PC is shown in Fig. 
|3la), which contains five periods of the index perturbation shown in Fig.|2f). The perturbation 
couples the incident waveguide mode upward, resulting in reduced waveguide transmission 
(< 0.5%) on the right side of the perturbation. In Fig. ISb), we consider the effect of this per- 
turbation on a higher-frequency mode at ft) = 0.3078, which also has even symmetry across the 
y = axis, but which has a periodicity of a in the direction of propagation. Equation (4) predicts 
that if Ae has a periodicity of 2a, the coupling constant fc^a will vanish in this case. This result 
is confirmed by a high in-plane transmission after the coupler (> 99%), as shown in Fig.[3jb). 
We will discuss this mode selective property of the coupler in greater detail in Section 3. 

Figure |4] displays the energy density of the y component of the evanescent field and the 
upward component of the scattering field from the coupler Figure |4ja) is a snapshot viewed in 
the vertical cross-section through a plane passing through the middle of waveguides, and Fig. 
Iltb) is the three dimensional view of the Ey energy density isosurface. We observe a strong 




I) 



ooooooooooo 

)OOOOOOOOOOC 

ooooooooooo 
OC 00 00 00 00 d 



DC DC DC DC DC 

oOoOoOoOoOo 

)OOOOOOOOOOC 

ooooooooooo 



15 5 



11 13 



10 
X 

(a) 



15 



20 



(c) 



Fig. 5. The control of the scattering direction, (a) Intensities of the scattered fields at dif- 
ferent feeding frequencies, (b) The structure of the a 5-period long coupler with increasing 
perturbation size. The brown shade area in the right panel labels the perturbation intro- 
duced, (c) The intensity distribution of the scattering field radiated from the coupler as 
shown in (b). The bottom panel is a x — z cross-section through the center of the first row 
of holes. 



radiated field originating from the 5-period long perturbed region. 

As observed in Fig. IH the radiation exits at 43° from the vertical. This direction depends 
on the interference of the fields scattered by discrete perturbations, as in a line grating. The 
interference is governed by relative amplitudes and phases of the fields at the positions of the 
perturbations. The scattered field can be approximated as a sum over all perturbations: 



i=l 

where 

^ ^(p^le'^O-f-^) (6) 

\r-n\ 

represents the propagation of the radiation from a scattering point at position ri, u{r) denotes 
the spatial profile of the radiated field, is the ^-vector in free space, ^v,, = {% /a+Niw,x)x is the 
crystal momentum of the perturbed PC waveguide where AA:„.,v indicates a possible shift in the 
A;- vector due to a slightly different dispersion in the perturbation region, ri is th position of the i- 
th scattering point, and n is the position of the start of the coupler The field dependence in Eq. 
(5) includes the decrease of the field amplitude along the coupler, a, and the phase difference of 
the waveguide mode between perturbations. Because the perturbations are n have a periodicity 
of lax, the phase difference is equivalent to laNin-^x- Equation (5) therefore predicts a shifting 



of the diffraction angle with wavelength, as given by the waveguide's dispersion relation. We 
observe this dependence in numerical simulations. In Fig. [3 a), the frequency of the incident 
waveguide mode is tuned from (0 = 0.2960 to CO = 0.3120 as the scattering angle shifts from 
-4° to 20°. 

The coupled mode theory discussed above assumes a weak perturbation to the waveguide 
mode. However, the mode intensity drops as the field propagates through the coupler. This 
decrease in mode intensity should be taken into account to improve the efficiency of the cou- 
pler Figure |3b) shows the structure of a 5-period long coupler in which we now include an 
increasing perturbation size to compensate for the decreasing field intensity. Using Eq. (5,6) 
to estimate the field intensity profile along the waveguide, we simulate a 5-period long cou- 
pler with perturbation sizes multipHed by factors of 0.8589, 0.8973, 0.9059, 0.9731, and 1. 
Figure |3c) shows the resulting intensity distribution of the scattered field viewed in a vertical 
plane through the center of the waveguides. We now observe a more uniform scattering along 
the length of the perturbation coupler, which results in higher mode overlap with the vertical 
Gaussian beam. 

3. Mode selectivity 

It is clear from the coupled mode picture that the out-coupling rate from the waveguide depends 
on the in-plane periodicity and the symmetry of the waveguide mode; for example, as was 
shown in Fig. [3] the in-plane transmission after the coupler is suppressed for the fundamental 
mode, but not for the higher-order mode at © = 0.3078. We now analyze this mode selectivity 
more closely in A:-space. Figure |6la) compares the dispersion relations for the fundamental 
mode (the blue line) and a higher mode (the green line) inside the photonic band gap, having 
frequencies of coa = 0.2436(27rc/fl) and cob = 0.3151 at = n/a. For the structure with index 
perturbations (the right panel), the dispersion is raised slightly because of the decrease in the 
effective refractive index. The periodicity of 2a scatters modes near A:, = tt/a by a lattice vector 
of Akx = n/ato the F point where = 0. Since these components are above the light line, they 
leak out of the PC plane and result in higher loss rate of the waveguide mode. 

In Fig.|6jb), the blue line shows the loss rate, y^, normalized by the loss rate of the fundamen- 
tal waveguide mode in the unperturbed structure, through different perturbation patterns. The 
green Une shows jr for the higher-order waveguide mode at o = 0.3151. As the perturbation 
size increases, the loss rate of the fundamental mode increases by over two orders of magni- 
tude, while the loss rate of the higher-mode remains roughly constant. This result demonstrates 
a high degree of selectivity of the perturbation coupler between different kinds of waveguide 
modes. Figure |6tc) shows the frequency increase of the fundamental and higher order mode 
(the blue and green lines). 

4. Transmittance and beam shape modification 

The efficiency of the coupler is estimated by calculating the transmittance from the waveguide 
mode to a target free space Gaussian spatial mode, which includes both upward and downward 
components. We approximate the efficiency by 

~rF^ \/p2\ ' ^ ' 

scatt I \^ Gl 

where S is the transmission out of the PC plane and the fraction quantifies the mode overlap 
of the upward(or downward)-traveling component of the radiated field with the target Gaussian 
beam profile. The scattering rate of the couplers with different index perturbations are shown 
in Fig. |7] (the red solid line). The perturbations not only scatter the incident waveguide mode 




Fig. 6. (a) Dispersion relations of waveguides without (left) and with index perturbation 
(right). The insets show the unit cells of the periodic structure for each case. The perturbed 
structure is the same as shown in Fig. Of), (b) Comparison of the changes in the relative 
decay rate of the fundamental mode (the blue line) and a higher mode (the green line), (c) 
Comparison of the changes in mode frequency of the fundamental mode (the blue line) and 
a higher mode (the green line). 



into the Gaussian mode, but also induce backward reflection (the gray solid line) and other 
in-plane leakage (the gray dotted line). The maximal scattering rate is reached when ro = 
and the perturbation size in the sides is 0.4466, in units of the fractional area of unit cell. As 
shown in Fig.|7] the scattering rate doubles when the perturbation size increases about 40% (the 
region A), but slightly decreases as the radius of the central holes increases (the region B). If a 
perturbation is introduced only at the center of the waveguide (the region C in Fig. |7]i, then the 
incident light is mostly reflected and only weakly scattered out of plane. 

When coupling to free space optics, and especially for single mode fiber coupling, the 
scattered field should have a large overlap with a Gaussian mode. In the case with the max- 
imal scattering rate of 0.6989, we calculate an overlap of 89% with a Gaussian beam fo- 
cused in the plane of the PC. This indicates that the coupler allows for both strong scat- 
tering and a highly Gaussian-like profile. The resulting transmittance through the coupler is 
Tr = 0.6989 X 89% = 62%. 

It is possible to decrease the beam divergence to allow for coupling to lower-NA lenses. By 



increasing the extent of the PPCWC in both x and y (transverse in-plane direction), we expand 
the coupHng region and decrease the output beam divergence. This modified index perturbation 
is shown in Fig.[8la), again optimized for the mode at to = 0.2436. Figures[8jb) and[8lc) show 
the intensity distributions of beams scattered from a small coupler (Fig. [3]) and a broadened 
one (Fig.[8ta)), viewed in y-z plane through the centers of the couplers. Figures |8ld) and[8te) 
compare the divergence of the radiated beams in both cases by fitting a Gaussian beam profile 
to the data. We assume 29 = 2Wo/zo, where 29 is the divergence angle of the beam and Wo and 
Zo are the beam waist and depth of the focus, respectively. In the case of the broadened coupler, 
the divergence is decreased by 20% from that of the original smaller coupler. 

Unfortunately, the out-of-plane scattering from the PPCWC is symmetrically upward and 
downward, so both contribute to the total scattering rate. While large Tc can be achieved by the 
methods described above, the overall efficiency of the PPCWC is so far limited to 50%. The 
upward scattering efficiency is defined as the directionality. 



D 



Pu 



Pu+Pd 



(8) 



where Pu and Po are the powers radiated in the upward and downward directions, respectively. 
One method of increasing D is to redirect the downward scattered field with a mirror This 
may be achieved with a DBR structure, separated from the PC slab by some distance, d. When 
d = 0.4667/;, the downward radiated light reflected by the mirror interferes with the upward 
radiated light constructively and a local maximum of D = 70% occurs. When we move the 
mirror further from the slab [d > 0.6/1), there is minimal interference between two beams due 
to divergence of the beams. The upward scattering rate is close to the total scattering rate in this 
case, and we achieve D = 95%. More sophisticated methods may be used to increase directional 
scattering, such as only partial etching of the PPCWC. This has been successfully employed 
for photonic crystal cavities ifTSlfTTl and channel waveguides 



5. Conclusions 

We have designed and implemented a method for selectively coupling a planar photonic crystal 
waveguide mode to a target free space mode. The method proceeds as follows: (1) one chooses 
a waveguide mode to out-couple, (2) one chooses an index perturbation to couple the mode that 
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Fig. 7. The calculated relative strength of transmission (black solid line), scattering (red 
solid line), reflection (gray solid line), and planar leakage (gray dotted line) of the 5-period 
long couplers with different index perturbation distributions introduced. The horizontal axis 
denotes the perturbation size, ro is the radius of holes at the center of the waveguide, and 
the perturbation size introduced on the side is in units of the fractional area of unit cell. 
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Fig. 8. (a) The structure of a coupler with additional index perturbation introduction to 
the 6th row from the waveguide, (b), (d) The intensity distribution viewed in a transverse 
plane through the center of the coupler shown in Fig. [3] (c), (e) The intensity distribution 
of the scattering from a coupler as shown in (a). A is the beam width (FWHM) measured at 
different heights above the center of waveguide, in units of lattice constant, a. 



has the appropriate periodicity and symmetry given by coupled mode theory, and finally, (3) one 
optimizes the perturbation for the desired out-coupled beam parameters. To demonstrate this 
approach, we focused on the problem of coupling a Wl waveguide mode to a slowly diverging 
Gaussian mode. We were able to ensure that the target waveguide mode was coupled with high 
rate, selectivity, and overlap with the target Gaussian mode. 

This perturbative photonic crystal waveguide coupler (PPCWC) not only enables efficient 
scattering of selected waveguide modes, but also allows modifications in the shape and the 
direction of scattering. The angle of the scattered field can be shifted by detuning the mode 
frequency with respect to the lattice constant or introducing some modulation to the perturba- 
tion size. By breaking the vertical symmetry of the PC or using reflectors under the PC slab, 
it is possible to minimize downward propagation, directing most vertical scattering upward. 
The PPCWC enables a straightforward integration of vertical in/out- couplers which can be 
employed in a range of application from optical interconnects to photonic crystal light sources. 
Furthermore, the spectrally resolved light coupler could find application in light trapping for 
thin-fibn solar photovoltaics. 
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